First order Mott transition at zero temperature in two dimensions: 
Variational plaquette study 
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PACS 71 . 10 . -w - Theories and models of many-electron systems 

PACS 71.30.+h - Metal- insulator transitions and other electronic transitions 

Abstract. - The nature of the metal- insulator Mott transition at zero temperature has been 
discussed for a number of years. Whether it occurs through a quantum critical point or through a 
first order transition is expected to profoundly influence the nature of the finite temperature phase 
diagram. In this paper, we study the zero temperature Mott transition in the two-dimensional 
Hubbard model on the square lattice with the variational cluster approximation. This takes into 
account the influence of antiferromagnetic short-range correlations. By contrast to single-site 
dynamical mean-field theory, the transition turns out to be first order even at zero temperature. 
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Introduction. — The correlation-driven transition 
from a paramagnetic normal Fermi liquid at weak cou- 
pling to a paramagnetic Mott insulator at strong coupling 
is one of the most important paradigms in solid-state the- 
ory [1,2]. For example, the Mott state is suggested to 
represent the proper starting point for theoretical stud- 
ies of the extremely rich and difficult correlation physics 
of two-dimensional systems such as cuprate-based high- 
temperature superconductors [3]. A big step forward in 
the understanding of the Mott transition was made by 
applying the dynamical mean-field theory (DMFT) [4-6] 
to the single-band Hubbard model which is believed to 
capture the main physics of the Mott transition in a pro- 
totypical way. The Hamiltonian is given by 



H = -* E c *W + u 



^2 ni t n a 



where the summation is over site indices 



the spin label, cj a and Cj i<r are the particle creation and 
annihilation operators and m t a = c\ a a a . Each doubly 
occupied site costs an energy U. All the numerical results 
are presented in units where t = 1. 

A number of objections concerning the DMFT picture of 
the Mott transition have been raised [7] and answered [8] . 
However, the main criticism comes from the neglect of 
the feedback of non-local magnetic correlations on the 
single-particle dynamics. This leads to a description of 
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and a is 



the paramagnetic Mott insulator with a macroscopically 
large ground-state entropy of S = L In 2 (where L — > oo is 
the number of sites). On the level of one-particle excita- 
tions, superexchange does not lift the 2 L -fold degeneracy 
arising from the spins of the localized electrons. This must 
be considered as a mean-field artifact which has profound 
consequences for the phase diagram: At finite tempera- 
tures, T > 0, the large entropy term stabilizes the Mott 
insulator in a situation where it competes with a metallic 
Fermi liquid. 

Cluster extensions of the DMFT [9] can cure this defect 
since they incorporate the feedback of short-range anti- 
ferromagnetic correlations. This has motivated a num- 
ber of previous studies which put the relevance of the 
mean-field picture for the transition into question and fo- 
cus on the qualitative change of the phase diagram due 
to short-range magnetic correlations in, say, two dimen- 
sions [10-12,14-17]. 

Recently, novel quantum Monte-Carlo techniques to 
solve the problem for a plaquette of correlated sites have 
been used to study the transition at finite T on the square 
lattice with nearest-neighbor hopping t. While different 
embeddings of the plaquette were considered, namely the 
cellular DMFT (C-DMFT) [16] and the dynamical cluster 
approximation (DCA) [17], it turned out that salient fea- 
tures of the DMFT phase diagram [5] are preserved (see 
fig. 1 ) . In particular, there is again coexistence of a metal- 
lic and an insulating phase in a certain U — T range at 
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Fig. 1: Top left: sketch of the DMFT phase diagram. Coex- 
istence of a metallic and a Mott insulating solution is found 
in the yellow region. A first-order Mott transition takes place 
at U C {T) (red line). The first-order line ends in second-order 
critical points (red dots) at T = T c and T = 0. Top right: pos- 
sible phase diagram within plaquette DMFT. The first-order 
line ends in U c i for T = 0. Bottom: the phase diagram sup- 
ported by our plaquette VCA calculations. The first-order line 
does not end in a critical point at T = 0. 

half-filling bounded by lines U c i(T) and U C 2(T). Compar- 
ison of their respective free energies leads to a first-order 
transition line U C (T) which, at a temperature T c , ends in a 
second-order critical point above which there is a smooth 
crossover only. 

For T < T c there are different possibilities: Within 
single-site DMFT (fig. 1, top left), the insulator wins at 
higher temperatures due to the residual high entropy of 
the insulator. This implies that U C (T) is increasing with 
decreasing T. Consequently, the line ends for T = in an- 
other second-order critical point which must coincide with 
the point up to which a metallic solution can be found, i.e. 
U c (0) — U C 2(0). At T — the metal is stable in the entire 
coexistence region. Contrary, within a plaquette DMFT, 
the entropy is low in the insulator due to short-range sin- 
glet formation. The phase diagram fig. 1 (top right) is 
obtained if the metal always wins at higher T. As com- 
pared with single-site DMFT, the critical line U C {T) bends 
back, and U c decreases with decreasing T [16]. Hence, 
U c (0) = E/ci(0), and at T = the insulator is stable in the 
entire coexistence region. 

Main results. — Here we use the variational cluster 
approximation (VCA) [18, 19] to embed a plaquette of four 
correlated sites and four or eight uncorrelated bath sites in 
the lattice. As the method is thermodynamically consis- 
tent and focuses on the optimization of a thermodynamical 
potential, it is ideally suited to distinguish between differ- 
ent phase diagram topologies. Here we consider T = 
using the Lanczos method as a cluster solver. Our re- 
sults are consistent with the previous plaquette DMFT 
studies [16, 17] but support yet another low-temperature 




Fig. 2: Reference systems consisting of four correlated sites 
with U ^= (filled, blue) and four (A) or eight (B) uncorrelated 
sites with U — (open, red). One-particle variational param- 
eters: "hybridization" V, hopping between correlated (t') and 
uncorrelated sites (t") and shift (±Ae) of the energies of the 
bath sites with respect to ^i. Here, fj, is the common chemical 
potential for both cluster and bath. For arbitrary one-particle 
parameters a space of trial self-energies £ is spanned on which 
a stationary point of the SFT grand potential Q[S] is searched. 

phase diagram (fig. 1, bottom). We find a sizable in- 
teraction range U c \ < U < U C 2 where the metallic and 
the insulating solution are coexisting at zero temperature. 
The T = endpoint of the line of first-order transitions 
U c = U C (T = 0) does neither coincide with U c \ nor with 
U C 2; so that the Mott transition is discontinuous also at 
T = 0. 

Method. — Using a plaquette of four sites to generate 
an approximate self-energy T,ij(u>) for the infinite D = 2 
square lattice, represents the essential step to go beyond 
the single-site DMFT approximation Sij(w) ~ 
An in principle ideal embedding of the cluster in the in- 
finite lattice could only be achieved with a continuum of 
bath degrees of freedom (uncorrelated sites with U = 0) 
attached to the correlated four-site cluster in the spirit of 
quantum-cluster theories [9]. For T = 0, however, this is 
not yet accessible with presently known cluster solvers. As 
far as static quantities and the thermodynamical phase di- 
agram are concerned, however, it is fortunately well known 
that a few bath sites can be sufficient for reliable predic- 
tions [20-23] . This holds on the single-site level where the 
DMFT phase diagram for the Mott transition can be re- 
covered qualitatively with a single bath site only [21] as 
well as for cluster approximations, as has been demon- 
strated recently in one dimension for the filling-controlled 
transition [22, 23]. 

To address the phase diagram of the Mott transition, we 
have to employ a thermodynamically consistent method 
to optimize the cluster self-energy and to fix the plaque- 
tte and bath parameters. For the reference systems dis- 
played in fig. 2 this can be achieved within the self-energy- 
functional theory (SFT). The SFT has been described in 
detail in refs. [18, 19, 21]. The main idea is to use the ref- 
erence system for spanning a space of trial self-energies. 
Referring to the reference system in fig. 2 (A), for example, 
the trial self-energy is parametrized as S = Sy,t',t" and 
calculated at a given U and at T = for each set (V, t' , t") 
from the cluster Green's function G' via Dyson's equation. 
On this space, the general SFT grand potential as a func- 
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Fig. 3: SFT grand potential shifted by the constant /iN, 
£"o[S] = Q[S] + fiN, and evaluated for trial self-energies Yly,o,o 
(see fig. 2 (A)) as a function of V. Results for different U at 
T — and for half-filling N = L (chemical potential fi = U/2). 
Circles indicate stationary points. Here, the value of the func- 
tional equals the ground-state energy Eo (per site). Filled cir- 
cles: stable metallic (blue) and insulating (red) phase. Shaded 
circles: metastable phases. At U c (arrow) there is a discontin- 
uous metal-insulator transition. A third (metastable) metallic 
phase (green circles) continuously coalesces with the insulating 
one at a hidden (metastable) critical point (hep). The energy 
unit is fixed by setting t = 1. 



tional of the self-energy, f2[S], can be evaluated exactly 



tt[£] = Q' + Tr In G - Tr In G' 



(2) 



Here ft' is the grand potential of the reference system 
and G is the lattice Green's function obtained via G = 
1/(Gq 1 — S) from the free lattice Green's function Go- A 
physical state is found as a stationary point of f2[Sy,t',t"] 
as a function of the variational parameters (V,t',t"). To 
this end we exploit the particle-hole symmetry of the Hub- 
bard model (1) at half- filling, use the (band) Lanczos 
method for a simultaneous calculation of all elements of 
the cluster Green's function G', and employ the Q-matrix 
technique or integration along the imaginary frequency 
axis to evaluate Tr In G/G' [24]. In addition, we make use 
of spatial symmetries to limit the number of independent 
one-particle parameters to a minimum. 

Parameter optimization. — For a simultaneous in- 
dependent optimization of V, if and t" (reference system 
fig. 2 (A)) we apply the downhill simplex method and/or 
iterative one-dimensional optimizations to find local min- 
ima of |Vf2[IV.t',t"]| 2 from which (if there are more than 
one) only those are retained for which il[S\/.t',t"] has a 



vanishing gradient. Particle- hole symmetry fixes the op- 
timal values e Ci0 pt = and £f,, op t = A* = U/2 for the 
on-site energies of the correlated and the bath sites, re- 
spectively. It has also been checked numerically that the 
SFT functional is stationary at these values. For the op- 
timal hopping between bath sites we find |i" pt | < t/25 in 
the entire U range considered. Note that for a finite num- 
ber of bath sites Lb the inclusion of t" enlarges the space 
of trial self-energies while for continuous baths (Lb — > oo) 
the approximation becomes equivalent [19] with C-DMFT 
where a coupling of baths attached to different correlated 
sites is not needed for this lattice geometry. 

The optimal hopping parameter between the correlated 
sites turns out as t' t = t + At' t with a small positive 



At' opt < t/10 for the U range considered here. A consid- 
erably larger At' opt is only found in the limit U — ► 0. For 
more itinerant electrons a stronger enhancement of the 
miracluster hopping is needed to (partially) compensate 
for switching off the mtercluster hopping in the approx- 
imation for the VCA self-energy. Note that t' opt = t for 
Lb — > oo. This is easily derived by a 1/w expansion of 
the SFT Eulcr equation (equivalently the C-DMFT self- 
consistency equation). The optimization of t' and t" ac- 
tually turns out to be almost irrelevant as compared to 
V. Setting t' = t and t" = and performing a one- 
dimensional optimization of V only, leads to changes in 
the optimal V of less than 1%. The change in the ground- 
state energy is negligible. 

Much more important is the inclusion of additional bath 
sites. While the reference system fig. 2 (A) with a single 
bath site (per correlated site) at the Fermi edge is expected 
to favor the metallic state, a slight bias towards the insu- 
lator is given with reference system (B) where two bath 
sites are taken into account at energies shifted by ±Ae 
away from the chemical potential /i = U/2. In case of 
(B), Ae and V are considered as independent variational 
parameters. 

Results and discussion. — We first concentrate on 
reference system (A) . Fig. 3 shows the dependence of the 
SFT functional on V. For U < U c2 « 6.35 (in units of 
t = 1) we find a metallic phase with a comparatively large 
optimal hybridization V op t which decreases with increas- 
ing U. For U > U c i « 4.6 there is a stationary point of the 
functional with a much lower V^ p t which is less U depen- 
dent. This corresponds to the Mott insulating phase as is 
obvious from the local Green's function and the self-energy 
displayed in fig. 4 for an interaction U = 5.8 in the coexis- 
tence region U c i < U < U C 2'. After Fourier transformation 
of the self-energy on the cluster, we find the self-energy 
for K = (0, 0) to be regular while for the cluster momen- 
tum K = (n, 0) it develops a pole at u = (fig. 4, right). 
This leads to a vanishing local Green's function Gn(iuj) 
for uj — > in the insulating phase (fig. 4, left). The same 
qualitative behavior has been seen at finite T in a recent 
C-DMFT study [16]. The metallic solution at U = 5.8 
is characterized by a finite Gn(iuj) and a regular Eif(iw) 
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Fig. 4: Imaginary part of the local Green's function Gu (left) 
and of the self-energy Ek (right) at the cluster momenta 
K — (0, 0) and K = (n, 0) as functions of imaginary frequen- 
cies iw at U = 5.8, i.e. in the coexistence regime, for the metal- 
lic (blue), the insulating (red) and the third phase (green). 
Results are obtained using reference system (A). Note that 
ImE( Tj o) = ImS(o,x) an d ImE(^ j7r ) = ImE( 0j o) due to particle- 
hole symmetry. 

for a; — > (fig. 4). Note that the K dependence is much 
weaker in the metallic solution. 

Comparing the ground-state energies at the respective 
stationary points (fig. 3), we find a discontinuous metal- 
insulator transition at a critical value U c = 5.79. The 
same picture and almost the same value for U c (within 
less than 0.1%) is found for the simultaneous and inde- 
pendent three-parameter (V, £',£") optimization. Also a 
different tiling of the square lattice (still using four-site 
plaquettes) or using a cluster with two correlated and two 
bath sites only (thereby breaking rotational symmetry) 
does not yield a qualitatively different picture. 

The critical U c from our plaquette VCA is substan- 
tially smaller than from single-site DMFT. The values 
t/i DMFT > « 11 (rcf. [12]) and C/ C (DMFT) « 12 (rcf. [17]) 

compare well with our mean-field result £/j DIA ' ~ 11.3 
which is obtained within the SFT by embedding a single 
correlated site coupled to a single bath site into the square 
lattice. Note that this two-site dynamical impurity ap- 
proximation (DIA) is known [21] to reproduce the DMFT 
phase diagram topology (fig. 1, top left). We have also 
verified numerically that it exhibits the entropy problem. 
Contrary, the plaquette VCA yields a vanishing ground- 
state entropy for (the metal and for) the Mott insula- 
tor and consequently supports quite a different picture, 
namely (fig. 1, bottom), a first-order transition at T = 0. 

Fig. 3 also demonstrates the presence of a third station- 
ary point (green) in the coexistence region which smoothly 
links the metal to the insulator but represents a mctastablc 
phase because of its higher ground-state energy. Note that 
in VCA, all stationary points are acceptable solutions. As 
can be verified from the finite quasiparticle weight and the 
finite density of states at the Fermi edge, this third phase 



is metallic in the entire coexistence region. This implies 
that for four bath sites, there is a quantum critical point 
at U c i which marks a continuous metal-insulator transi- 
tion. The additional stationary point is reminiscent of the 
third solution which is found [25] in single-site DMFT at 
finite temperatures and which also interpolates between 
the two main phases. At finite T, however, this is more 
like a gradual crossover opposed to the (hidden) quantum 
critical point obtained here. 

Formally, the plaquette VCA for a single-band model 
becomes equivalent with the single-site DIA applied to a 
model with four orbitals per site, namely if the four cluster 
momenta are identified with the four orbitals. Contrary 
to the doped system [26] , the continuous transition at U c \ 
is not orbital selective (in this interpretation), i.e. with 
decreasing U all four A'-points simultaneously undergo 
the transition to the insulator (develop a gap in the in- 
dependent spectral function). The complicated structure 
of the Coulomb interaction in K space, and the presence 
of corrclated-hopping (inter-orbital) terms in particular, 
makes orbital selectivity implausible. For K = (it, 0) and 
K = (0, 7r) the transition is of the Mott-Hubbard type. 
This is consistent with the fact that for the particle-hole 
symmetric case and in the Mott-insulating state, the self- 
energy at lu = diverges on the non-interacting Fermi 
surface [27]. For K = (0,0) and K = (n,Tr), the insu- 
lating spectral function is neither Mott-Hubbard-like (as 
it is not particle-hole symmetric) nor band-insulator-likc 
(as the "orbital" K = (0, 0) is not fully occupied and the 
"orbital" K = (tt,tt) not completely empty). In the con- 
text of multi-orbital DMFT, a transition to an insulating 
state with almost complete orbital polarization has been 
discussed for the titanates [28]. 

We now move to the reference system fig. 2 (B) with 
two bath sites per correlated site. Again one finds a first 
order transition, but with different critical values for the 
interaction strengths. To be more specific, fig. 5 shows 
the ground-state energy and the double occupancy. There 
are two sets of disconnected solutions. Since the metal- 
lic and insulating solutions with respective lowest energy 
belong to families that are not directly connected, it be- 
comes meaningless to try to find the analog of the hidden 
critical point obtained for reference system (A). 

As compared to (A), the coexistence range shrinks. We 
define the coexistence range from hysteresis. In other 
words, decreasing U starting from large U, we follow the 
insulating solution (red line) until it disappears. This de- 
fines U c \ = 5.25. We then follow the metallic solution 
(blue line), increasing U from small U, until it disappears 
at U C 2 = 6.37. These two numbers are more clearly seen 
on the double-occupancy plot on the bottom panel of fig. 
5. The difference AU coex . = U C 2 — U c \ = 1.12 is not very 
different from the T = extrapolation of the finite-T C- 
DMFT results of ref. [16] (A[/ cocx . « 0.73) and clearly 
improves the result AU cocx , = 1.75 obtained with (A). 
This suggests convergence with respect to the number of 
bath sites similar to the single-site DIA [29]. Note, how- 
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Fig. 5: U dependence of the ground-state energy Eq (top 
panel) and the double occupancy {rufriii) (bottom panel) in 
the metallic (red) and in the insulating (blue) phase with refer- 
ence system (B). The two solutions are not connected; however, 
the insulating solution is smoothly connected to an unstable, 
and then metallic solution (dashed curves) that together form 
a pattern similar to the one observed with reference system 
(A). That metallic solution has however a higher energy than 
the disconnected metal solution (blue). The latter is also con- 
nected to an unstable solution (not shown) very close to it. The 
coexistence range, U c i and U C 2, is indicated. The actual tran- 
sition takes places at U c where the energies of the two solutions 
cross. 



ever, that our results differ from those of C-DMFT solved 
by exact diagonalization at T = [f 2, 13]. 

The critical U c for the actual metal-insulator transition 
can be read off from fig. 5 as U c — 5.5 to be compared 
with U c = 5.8 obtained with (A). As could be expected 
from previous cluster studies [22], the reference system (B) 
favors the insulating phase and brings U c closer to U c i as 
compared to (A). A reference system with three bath sites 
would favor the metallic phase again. We therefore believe 
that the transition would remain first order even if more 
bath sites were added. 

Within DMFT there is an insulating solution above and 



below U C 2 from which at U C 2 a metallic solution splits off 
when decreasing U. This is a bifurcation of the insulating 
solution of the non-linear DMFT equation. By contrast, 
within the plaquette approach such a bifurcation mech- 
anism is not necessary since the metal and the insulator 
have the same (vanishing) ground-state entropy. 

Discussion and conclusion. — The Mott insulator is 
best characterized at high temperature by two properties: 
a) insulating behavior, i.e., a gap or pseudogap between 
two peaks in the single-particle density of states, b) no 
long-range order [30] . In the Mott insulator at fixed high 
T, if one decreases U, one crosses over to a metallic-like 
state where the density of states has a single peak. One 
expects the crossover to be replaced by a phase transition 
at low temperature. In reality, long-range antiferromag- 
netic order is present in the ground state of the Hubbard 
model with nearest-neighbor hopping only. At infinitesi- 
mal U it is driven by nesting. In other words, at small 
U the gap originates from Slater physics where antifer- 
romagnetic correlations increase with U and the gap in- 
creases as exp (—2iry/ (t/U)). At very large U, one is in 
the Heisenberg limit where the antiferromagnetic correla- 
tions decrease with U and the gap increases linearly with 
U, so Mott physics is relevant. And there is no phase 
transition between these two limiting antiferromagnets. A 
phase transition occurs at low enough temperature if an- 
tiferromagnetic long-range order is prohibited. In single- 
site DMFT, one cuts off the correlation length to zero and 
there is a T — second-order transition at U about 1.5 
times the bandwidth. If one lets the antiferromagnetic 
correlations grow, the longer the range of these correla- 
tions, the closer the transition will be to U = 0. In the 
case we studied, antiferromagnetic correlations are cutoff 
beyond second neighbor. We suggest that this case is the 
one that is most closely connected to the high-temperature 
crossover that characterizes a Mott insulator. Indeed, the 
value U c = 5.5 that we find for the first order transition is 
closest to the one extrapolated from the finite temperature 
crossover studied in quantum Monte-Carlo calculations on 
finite lattices with up to 8 x 8 sites [31]. 

At finite T, first order transitions are best understood as 
a tradeoff between energy and entropy in different phases. 
At T = 0, we should consider the tradoff between poten- 
tial and kinetic energy in various phases as the mecha- 
nism for the transition. An insulating phase, whether it 
is stable or metastable, always has lower potential energy 
than a metallic phase. However, between U c \ and U c the 
lower kinetic energy (less localization) of the metal makes 
it more stable than the metastable insulator despite the 
potential energy advantage of the insulator. Between U c 
and U C 2, the kinetic energy of the metastable metal is 
not small enough compared with that of the insulator to 
overcome its lower potential energy. Indeed, the kinetic 
energy of the insulator can at best grow like —t 2 /U in the 
large U limit while the kinetic energy of the metal can 
grow faster as electrons are scattered further away from 



Matthias Balzcr, 1 Bumsoo Kyung, 2 David Senechal, 2 A.-M. S. Tremblay, 2 Michael Potthoff 3 



the Fermi surface by the interaction. In practice, the first 
order transition that we found is in the intermediate cou- 
pling regime where the T — insulator is neither clearly 
in the Heiscnbcrg (Mott) limit nor in the Slater limit. 

The central idea of dynamical (cluster) mean-field the- 
ory in general is that the analysis of the quantum critical 
point in the paramagnetic state provides the key to an 
understanding of the entire phase diagram - even if this 
point is obscured by long-range magnetic order: Namely, 
aT = critical point of a non-magnetic origin will also 
rule the physics above the ordering temperature, in the 
doped system or in the presence of magnetic frustration. 
In the past, this concept has been frequently used at the 
single-site DMFT level. 

Our main result, however, is that such a quantum criti- 
cal point is absent. The T = transition is first order. The 
absence of a continuous non-magnetic T = transition is 
contrary to the widespread DMFT result and therefore ex- 
pected to have profound consequences for our understand- 
ing of (doped) Mott insulators in two dimensions. It has 
been argued based on Hartree-Fock calculations that ex- 
plicitly introducing frustration through longer-range hop- 
ping should restore quantum critical behavior at finite 
frustration [32], but that remains to be verified. 
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